Magneto- Acoustic Waves of Small Amplitude 
in Optically Thin Quasi-Isentropic Plasmas 



Valery M. Nakariakov 

School of Mathematical and Computational Sciences, University of St Andrews, St Andrews, Fife KY16 

9SS, Scotland, UK 

and 

Cesar A. Mendoza-Briceno, Miguel H. Ibanez S. 

Centro de Astrofisica Teorica, Universidad de los Andes, Apartado de Correos 26, IPOSTEL La 

Hechicera, Merida, Venezuela 



ABSTRACT 

The evolution of quasi-isentropic magnetohydrodynamic waves of small but finite 
amplitude in an optically thin plasma is analyzed. The plasma is assumed to be initially 
homogeneous, in thermal equilibrium and with a straight and homogeneous magnetic 
field frozen in. Depending on the particular form of the heating/cooling function, the 
plasma may act as a dissipative or active medium for magnetoacoustic waves, while 
Alfven waves are not directly affected. An evolutionary equation for fast and slow 
magnetoacoustic waves in the single wave limit, has been derived and solved, allowing 
us to analyse the wave modification by competition of weakly nonlinear and quasi- 
isentropic effects. It was shown that the sign of the quasi-isentropic term determines 
the scenario of the evolution, either dissipative or active. In the dissipative case, when 
the plasma is first order isentropically stable the magnetoacoustic waves are damped 
and the time for shock wave formation is delayed. However, in the active case when the 
plasma is isentropically overstable, the wave amplitude grows, the strength of the shock 
increases and the breaking time decreases. The magnitude of the above effects depends 
upon the angle between the wave vector and the magnetic field. For hot (T > 10 
K) atomic plasmas with solar abundances either in the interstellar medium or in the 
solar atmosphere, as well as for the cold (T < 10 3 K) ISM molecular gas, the range 
of temperature where the plasma is isentropically unstable and the corresponding time 
and length-scale for wave breaking have been found. 



Subject headings: Magnetoacoustic Waves, Thermal Instability, Shock Waves, Interstel- 
lar Medium, Solar Atmosphere 
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1. Introduction 

The study of wave propagation and the sta- 
bility of optically thin plasmas is of crucial im- 
portance for understanding the origin and evolu- 
tion of inhomogeneities observed at very differ- 
ent length scales and time scales in astrophysical 
plasmas. Travelling acoustic waves in the plasma 
can be unstable due to peculiarities of the pro- 
cesses of heating and cooling of the medium. The 
possible amplification of the waves is connected 
with additional energy liberation in the plasma 
in the compression phase or, in other words, with 
nonadiabaticity of the plasma. In different astro- 
physical applications, the nonadiabaticity of the 
plasma can be mathematically modelled by an 
introduction of an additional term in the adia- 
batic equation. This term, usually called a heat- 
ing/cooling function, depends on parameters of 
the plasma under consideration. The specific ex- 
pression for this function varies in different phys- 
ical situation (Rosner, et. al., 1978; Vesecky et. 
al., 1979; Dahlburg and Mariska, 1988), but for- 
tunately some general properties of this function 
define main features of dynamics of perturbations 
in optically thin plasmas (either amplification or 
decay) . 

The above problems, taking into account many 
different physical effects have been extensively 
studied as far as the linear approximation is con- 
cerned (Parker, 1953; Kruskal and Schwarzschild, 
1966; Weymann, 1960; Field, 1965; Priest, 1982; 
Ibahez and Parravano, 1983; Ibahez, 1985; Ibaahez 
and Escalona, 1993; Bodo et al.,1985; Balbus and 
Soker, 1989 and references therein). The linear 
approach corresponds to the case of waves with 
zero amplitude and, consequently, can describe 
only an initial stage of the instabilities associated 
with nonadiabatic effects in the plasma. Gener- 
ally, the nonlinear regime is not far well under- 
stood and most of the progress has been done by 
numerical simulations of very particular struc- 
tures (Dahlburg et al., 1987; Karpen et al., 1989; 
Reale et al., 1994). 



In addition to condensations that can be formed 
in hot plasmas by thermal instability in a nonadi- 
abatic plasma, providing that the instability does 
not saturate in the nonlinear regime, the steeping 
of travelling waves may also originate inhomo- 
geneities (Oppenheimer, 1977; Krasnobaev, 1975; 
Krasnobaev and Tarev, 1987 (Reference KT herein 
after); Tarev, 1993). This problem has been 
worked out in KT and Tarev (1993) references 
neglecting the presence of magnetic fields. How- 
ever, in many astrophysical situations the mag- 
netic fields become relevant, therefore, it is wor- 
thy to analyse the propagation of nonlinear mag- 
netoacoustic waves in plasmas where nonadia- 
batic effects also go into play. This problem be- 
yond its heuristic interest, also can be the clue 
for understanding particular observed features in 
magnetised plasmas, for instance, the intensity 
oscillations observed in the network bright points 
in the quiet Sun (Kalkofen, 1997). Also, in the 
presence of the magnetic field, the magnetoacous- 
tic waves perturb the absolute value and com- 
ponents of the field, which gives an additional 
possibility for the registration of the discussed 
phenomena in astrophysical plasmas through ob- 
servation of variation of gyro-emission, Zeeman 
effect etc. 

The aim of the present paper is the general 
study of weakly nonlinear magnetoacoustic and 
Alfven wave dynamics in presence of nonadia- 
batic effects in an optically thin plasma. We con- 
sider a homogeneous plasma penetrated by the 
straight and homogeneous magnetic field. Ef- 
fects of self-gravitation, ionization and steady 
flows are neglected. The physical state is deter- 
mined by time and spatial independent variables. 
Additionally, the heating/cooling function (per 
unit volume and time) is assumed to be depen- 
dent of pressure and density, i.e. Q(p,p). Only 
quadratic nonlinear effects will be taken into ac- 
count. Assumption that nonlinear and nonadia- 
batic effects are weak, but not necessarily of the 
same order, allows us to apply the method of 
slowly varying amplitudes and derive an evolu- 
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tionary equation for fast and slow magnetoacous- 
tic waves. Note that, despite the nonlinearity 
considered in the paper is weak, it leads to the 
significant modification of the wave, up to for- 
mation of shocks. In the derivation of the evolu- 
tionary equation, a single wave approximation is 
applied. This means that interaction of waves of 
different types, as well as interaction of forward 
and backward waves, is neglected. The nonlin- 
ear interaction of different modes, for example, 
three-wave interaction, occurs when all interact- 
ing modes are existing in the system (and are, 
actually, over some certain amplitude threshold 
prescribed by dissipation). Also, efficiency of the 
interaction depends strongly on coherence of the 
interacting waves. Both these conditions are as- 
sumed to be not fulfilled: we consider a wave of 
one certain type only, while other modes are not 
excited. Also, our attention is restricted to evo- 
lution of wave pulses with a wide spectrum, and 
nonlinear resonant interactions are suppressed. 
Thus, the dominating nonlinear process is gener- 
ation of higher harmonics, leading to wave steep- 
ening and shock formation. Moreover, this pro- 
cess is more pronounced due to the absence of 
linear high frequency dispersion. Nonlinear dis- 
persion preventing the second harmonics genera- 
tion, may appear in the next order of nonlinearity 
and is out of scopes of the paper. 

Since the thermally unstable plasma is an ac- 
tive medium for the magnetoacoustic waves, the 
solution of the evolutionary equation shows the 
amplification of the initial magnetoacoustic per- 
turbation and speeding up of shock formation. 
The breaking time depends on the angle of the 
wave propagation with respect to the magnetic 
field and is different for slow and fast waves. The 
above results are quite general and potentially 
can be applied to different astrophysical as well 
as laboratory plasmas. However, in this paper, 
the results obtained will be addressed to study 
a hot (T > 10 4 K) atomic plasma (Solar atmo- 
sphere or hot Interstellar Medium) and a cold 
(T < 10 3 K) molecular gas (say for instance, the 



cold ISM gas). 

2. Governing equations 

In the model considered, dynamics of the 
plasma is described by the set of MHD equations, 



dB 

~dt 



= curl V x B, 



„ 1 

p— — = - Vp — — B x curl B, 

at 47r 



^ + div P V = 0, 

dp jpdp 

^"7^ = (7_1) Q{p ' pl 



divB = 0, 



(1) 
(2) 
(3) 
(4) 
(5) 



where Q(p,p) is a known function describing ef- 
fects of the nonadiabaticity and all other nota- 
tions are standard. 

The stationary magnetic field is in the xz— plane, 
Bo = -Bosnia xo + -Bocosazo , where Bq is the ab- 
solute value of the magnetic field, a is the angle 
between the magnetic field and z— axis, xo and 
zo are the unit vectors. We consider dynamics 
of waves propagating along the z— axis. Depen- 
dences upon x and y are ignored (d/dx = d/dy = 
0). 

Projecting Eqs. (1) — (4) on the axes, we have 



dB x 
dt 



d_ 

dz 



(V Z B X - V X B Z ), 



dBz 
dt 



= 0, 



dV x _ 1 dB x 
' dt 4vr 2 dz ' 



(6) 

(7) 
(8) 
(9) 
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dVy 
) - 

d t 



An z dz 



dVz 
7 d t 



dz An V x dz v dz 



i + ifz^^ 



(10) 



(11) 



(12) 



dp dp ip (dp „dp\ 

(13) 

We rewrite equations @-fll3|) gathering linear 
and nonlinear terms on left and right handside, 
respectively, 



dB d 
dB„ 



£ (B 0z V y ) = N 2 , 



dt dz 
dV x B 0z dB x 



Po- 



An dz 



Po 



dVy B Qz dBy _ N 



Po- 



dt An dz 
dV z , dp , B 0x dB x 



+ 7^ + 



dt dz An dz 



dp dp (dQ dQ 



dp dv z 

dt dz 



(14) 

(15) 
(16) 
(17) 
(18) 

N 7 , 
(19) 

(20) 



where Bq x = Bq sin a and Bq z = Bq cos a; N\, 
N 2 , N^-Ng are nonlinear terms. The derivatives 
of the heating/cooling function Q in equation 
( |l9| ) are calculated for unperturbed values of the 
density and pressure, p = po and p = po- Accord- 
ing to equation @, we choose the perturbation 
B z = 0. 



Restricting ourselves to keep quadratic non- 
linear terms only, we obtain the following expres- 
sions for the nonlinear right handsides, 



d_ 

dz 



d 



(V Z B X ) 



N 2 = --(V z B y ), 



dv x dv x 

-P^r ~ PoV z 



dt 



dVy 

'-ft 



PoV z 



dz 

My 
dz 



(21) 

(22) 
(23) 
(24) 



dV z 



dV 7 1 „ dB T 1 „ dB„ 



(25) 



at d P m dp dp dp 

N7 = - p m- poVz d-z + ^dI +1PoVz d-z' (26) 



N S = -§- z (pV z ). 



(27) 



In Eqs. ( pl| ) — (p7|), variables p, p, V x and B x 
may be expressed through the variable V z from 
the following relations: 



dp _ dV z 
dt Po dz 



~dl = - poCs -d7> 



(28) 



(29) 



d 2 V ( d 2 d 2 \ 

^ = -» ta^-^ K (30) 



d 2 r2 2 d 2 



B, 



Bq sin a 



d 2 V z 
dzdt 



(31) 
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Eqs. ( |l5l) and (|i~7| ) can be combine to the equa- 
tion for the Alfven wave: 

1 (dN 5 Bo z dN 2 \ 

Vaz Vy = 7A~dT + 17^7j ' (32) 

where we use the Alfven wave operator, 

and C Az = S 0z /(4^o) 1/2 - 

Eqs. (P^) - (PC|) can be combined to the equa- 
tion for magnetacoustic waves: 



V Az V s V z 



r 2 d'V z 



: dz 2 dt 2 



(7-1) 



dQ +c 2dQ 
dp s dp 



dN* 



d 2 v z 

dz 2 



dt 



Po I 

Box d 2 



dt 



^ (c 2 N 8 + -N, 
dz V Po 



4vr dzdt V <9i 



(33) 



where 



dt 2 



c 



s dz 2 ' 



C s = ylPofpo is the sound speed and Cax = 
5ox/(4vrp ) 1/2 . 

In the linear limit, the Alfven wave perturbs 
the variables B y and V y , the magnetosonic waves 
perturb the variables B x , V x , V z , p and p. As it 
follows from (|32|), the introduction of the heat- 
ing/cooling function Q(p,p) does not change the 
propagation of the Alfven waves, at least till the 
quadratic nonlinear terms. 

Equations (^) and ([33]) form the governing 
set of equations for consideration of quadrati- 
cally nonlinear MHD wave dynamics in an homo- 
geneous plasma with nonadiabaticity. Set (p2|), 
(|33|) contains information about independent dy- 
namics of the MHD waves (Alfven and fast and 
slow magnetoacoustic) as well as their nonlin- 
ear interaction. This set is an analogue of the 



similar set of equations derived in by Nakari- 
akov and Oraevsky (1995), Nakariakov et al. 
(1997a,b, 1998) for MHD waves of finite ampli- 
tude in smoothly inhomogeneous cold plasma. 

An analysis of the right handsides of equa- 
tions (^) and ( |33| ) shows that if the Alfvenic 
perturbations (i.e. B y and V y ) are initially ab- 
sent from the system, they can not be excited by 
the magnetoacoustic perturbations, because all 
terms on the right handside of ( ^2|) are products 
of Alfven and magnetoacoustic variables. This 
gives an advantage for an analytical description 
of the nonlinear dynamics of the magnetoacous- 
tic waves, allowing to take V y and B y zero always 
and everywhere. Strictly speaking, this assump- 
tion is non-physical, because Alfvenic perturba- 
tions are always presented in the realistic plasma 
as noise. Consequently, in the natural condi- 
tions there is a nonlinear interaction of Alfven 
and magnetoacoustic waves. Nevertheless, if the 
level of Alfvenic noise is initially sufficiently low, 
the assumption can be applied to the descrip- 
tion of magnetoacoustic waves in some certain 
initial stage of the interaction, when nonlinear 
terms including Alfven variables are much less 
than magnetoacoustic terms. 

3. Dispersion relations 

Supposing V z ~ exp(iu;i — ikz), we obtain the 
dispersion relation for the magnetosonic waves 
from Eq. @: 

{uj 2 - C\ cos 2 a k 2 ){oj 2 - C 2 k 2 ) - C 2 A sin 2 aoj 2 k 2 

+ -Ak 2 (u? - C\ cos 2 a k 2 ) = 0, (34) 
u> 

where Ca = Bo/i^irpo) 1 ^ 2 and 



(35) 



If the heating/cooling function is not taken 
into account, A = 0, we have the well-known 
dispersion relation for magnetosonic waves, 

{oj 2 -C a cos 2 a k 2 ){oj 2 -C 2 k 2 )-C A sin 2 a u?k 2 = 0, 

(36) 
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which describes dispersion for both fast and slow 
magnetosonic waves. Dispersion relation ( |36l) 
gives that the phase speed is 



r 2 

^fast, slow 



1 



(c 2 s + c\y 



4 cos 2 a C 2 S C\ 



(37) 

where the upper sign corresponds to the fast wave 
and the lower to the slow wave. Figure (|l]) shows 
the variation of this speed as a function of a for 
different values of plasma beta (/?). 

Consider the case of the parallel propagation, 
a = 0, 



(yj 2 - C 2 A k 2 )(uj 



C 2 k 2 + -Ak 2 ) 
u> 



0. (38) 



There are fast and slow waves. The thermal in- 
stability can be only on the slow mode, in this 
case. If Ca = 0, we obtain the same result as in 
KT reference. If C s — > (the coronal case), we 
have a> 3 = —iAk 2 . 

Consider transversal propagation, a = w/2, 



2 / 2 

UJ UJ 



(C 2 A + C 



^k 2 

UJ 



0, (39) 



there is not slow wave now, but the fast wave can 
be a subject to the thermal instability. 

In the hydrodynamical case, Ca = 0, the dis- 
persion equation becomes 



UJ 



Cik 1 + -Ak 2 

UJ 



0, (40) 
reproducing the dispersion relation of KT. 



4. Derivation of the evolutionary equa- 
tion 

Now we take that the influence of the heat- 
ing/cooling function on the wave propagation is 
small, in other words, we suppose that in equa- 
tion ( |33|) the third term on the left handside is 
much less than the other linear terms. Also, 
the quadratically nonlinear terms are taken into 
account. In the weakly nonlinear, weakly non- 
adiabatic case, the magnetosonic waves propa- 
gate with the phase speed about the speeds given 



by expressions (37). The weak nonlinearity and 
non-adiabaticity lead to slow evolution of the 
waves. (The term "slow" means that the typical 
times of the nonlinear and non-adiabatic modifi- 
cation are much greater than the wave period) 

Following one certain magnetoacoustic wave, 
we change to the running frame of reference, 



£ = Z-Ct, T = t, 



(41) 



where C is either Cfast for the fast wave or C s i ow 
for the slow wave. The variable r means the slow 
time, which corresponds to the slow evolution 
of the wave considered due to the weak nona- 
diabaticity and nonlinearity. The dependence of 
the wave amplitude upon this slow time r rep- 
resents the difference of the wave propagation 
in the presence of the weak nonadiabaticity and 
nonlinearity from the ideal and linear case. In 
this frame of reference, equation (^) is rewrit- 
ten as 



-2C{2C 2 -C 2 -C 2 A ) 



V z +A 



C 2 



C\ cos 2 a d 3 



Po <9£ 3 
Bn sin a 

+ 



drde 
(C 2 - C\ cos 2 a) 



C 



i 



v. 



CN 6 + (C 2 S N S + —N 7 ) 
Po 



C 



-Cm + ^ C ° Sa Ar 4 



(42) 



47T \ po 

where expressions for the nonlinear terms are fol- 
lowing: 



N- 



B Csma dV 2 



(43) 



1 C 2 -C;|cos 2 a 8£ ' 

N 4 = 0, (44) 

B 2 C 2 sin 2 a dV z 
4^ [C 2 - C 2 A cos 2 a) 2 " &p 

C z [ ' 



(45) 



N 7 



Na 



2_po m 

C z d£ 



(47) 



Only first order derivatives with respect to the 
slow time r are kept on the left hand side of equa- 



tion (42) and neglected in the non-adiabatic and 
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nonlinear terms. We used expressions (g8[) -(pj ) 
rewritten in the moving frame of reference ( [111) > 



P=-gV z , p = p — V z , V x 



c 2 -c 2 r 

cot q — — — 14 



C 2 



B, 



BqC sin a 
C 2 - C\ cos 2 a 



V z 



(48) 



where only fast dependence on time has been 
taken into account. 

Using (|43| ) — (^7|), we can rewrite equation ( f42|) 

as 



a 4 



(9 3 (9V 

■V 8 + M-^^ + e-5^^-57 L = 0, (49) 



5^ 

^e 3 



d£ 3 d£ 



where fj, and e are the coefficients, 



/i 



-A 



C 2 



C\ cos 2 a 



2C 2 (2C 2 - C 2 



(50) 



_ 3CjC 4 sin 2 a + C 2 ( 7 + 1)(C 2 - C\ cos 2 a) 2 

£ ~ 2C 2 {C 2 - C\ cos 2 a)(2C 2 - Cf - C\) ' 

(51) 

Integrating equation ([i"9|), and assuming that 
the arbitrary constants of the integration are 
zero, we obtain an evolutionary equation describ- 
ing the weakly nonlinear dynamics of magnetoa- 
coustic waves in the magnetised plasma in the 
presence of the heating/cooling function, 



or a£ 



0. 



(52) 



Equation ( |52|) is an analogue of the Burgers and 
Korteweg - de Vries equations (the difference is 
connected with taking into account of the effect 
of nonadiabaticity instead of dissipation or dis- 
persion leading to appearance of the terms with 
highest derivatives with respect to the spatial co- 
ordinate) , and is a generalisation of the equation 
obtained in Ref. KT for acoustic waves. 

We note that equation ( |52| ) describes the weakly 
nonlinear and nonadiabatic evolution of either 
the fast magnetoacoustic wave or slow magnetoa- 
coustic wave, but not their nonlinear coupling. 



The nonlinear interaction of the magnetoacoustic 
waves is neglected due to the single wave approx- 
imation applied above. In this approximation, it 
is supposed that there initially is only a wave 
of one certain type in the system, while another 
wave is absent. Nonlinear coupling of the MHD 
waves in the thermal unstable medium will be 
considered elsewhere. 

Consider dependence of the coefficients /x and 
e of evolutionary equation (|32|) on the problem 
parameters. Actually, there are two independent 
parameters in the problem considered, the an- 
gle of the wave propagation a and parameter f3 
defined ratio of the kinetic and magnetic 
pressure in the plasma, (3 = ^C 2 /C\. 

Figures @ and @ show these two coefficients 
as a function of a for different values of /?, for fast 
and slow waves, respectively. In Figure (^) can 
be seen for a given value of the angle both coeffi- 
cients increase when (3 decreases but the opposite 
occurs for slow waves (see Figure (||)). When 
(3 = 1.2 both coefficients (^C 2 |A| _1 and e) are 
independent of angle a. 

In the case Bq = (Ca = 0, (3 — > oo), we have 



C — Cfast — C s ; [i 



Expressions (^) coincide with the expressions 
obtained in Ref. KT. All coefficients are isotropic 
(independent on the angle of the wave propaga- 
tion). 

5. Analysis of wave dynamics 



Equation ( p2[ ) is a partial differential equation of 
the first order, and can be easily solved analyti- 
cally in an implicit form (KT), 

V z -fU + -V z [1 - exp(/xr)A exp(-/ir) = 0, 

(54) 

where the function F{x) is a profile of the wave 
in the initial time, V z (z,0) = F{z). 

For convenience of the following consideration, 



we rewrite solution (54) in the laboratory frame 
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of reference, 

V z -F \ z-Ct + -V z [l- exp(jit)]j exp{-fd) = 0. 

" (55) 
When the parameter of nonadiabaticity fi tends 
to zero, solution (|55| ) becomes the well-known 
simple wave solution (Landau and Lifshitz, 1987) 



V z - F(z -Ct- eV z t) 







(56) 



describing steepening and breaking of a nonlin- 
ear wave in the absence of dissipation, disper- 
sion and nonadiabaticity. Nonlinearity leads to 
generation of second harmonics, carrying out the 
upward energy shift in the wave spectrum. In 
a finite time, there appear infinite gradients in 
the wave shape, corresponding to a weak shock 
wave. If the parameter of nonlinearity e is also 
zero (this takes place in the case of a wave with 



infinitely small amplitude), solution (56) shows 
that the wave propagate without a change of the 
wave shape. Different signs of the parameter e 
define where, at the front or rear slope of the 
wave (or, more precisely, on the slope with posi- 
tive or negative sign of derivative of the function 
describing the wave profile), the shock formation 
does take place. 

To obtain the time and position of formation 
of the shock wave, or the time and coordinate of 
where the wave is breaking, it is convenient to use 
another implicit form of the solution of equation 
( |52|) which can be easily obtained from solution 



z-Ct+-V z (l 



e M ) + f(V z e^) = 0, (57) 



where the function / is defined implicitly by the 
initial shape of the wave, z = f(V z ). 

Following Ref. KT, we find that the breaking 
time is given by the expression 



^breaking 



--log 

I* 



1 + */ 

e 



(58) 



where / is the value of the derivative of the 
wave profile at the inflection point. In the adia- 
batic case (/z = 0), formula (|58|) reduces to the 



well-known formula for a simple wave (Landau 
and Lifshitz, 1987). In the unmagnetised case 
Bq = 0, formula ( |58[ ) reduces to the formula for 
acoustic wave breaking time. Since the parame- 
ter e is positive for all angles of propagation for 
both fast and slow waves, the shock appears at 
the slope with / < for all angles and both 
waves. The specific value of the breaking time 
depends on the type of the wave considered and 
the angle of propagation. 

Modification of the wave evolution, carried 
by finite nonadiabaticity, depends on the sign 
of the parameter /i. A positive \x acts similarly 
to dissipation, leading to wave decay. In such a 
case, equation (|58|) defines a maximum value of 
A* = Aimax = £■/!/' | beyond which the damping of 
the disturbance by non-isentropic processes oc- 
curs in time shorter than the time at which non- 
linear effects go into play, and shocks waves can 
not be formed. A negative [i corresponds to the 
case of thermal instability, when the wave ampli- 
tude grows in time. (More strictly, in this case, 
this effect should be called thermal over- 

stability or oscillatory instability). According to 
( |50D and Fig. 2a and 3a, the sign of the coefficient 
of nonadiabaticity is opposite to the sign of the 
parameter A introduced in (|35"|). Consequently, 
the amplification of the finite amplitude magne- 
toacoustic wave takes place when the parameter 
A is positive. 

Different regimes of evolution of the initial 
Gaussian pulse, 



V z (z) = V exp(-z 2 /L 2 ), 



(59) 



where Vq is an amplitude and L is a characteris- 
tic length of the pulse in the initial time t = 0, 
are shown on Figure 4. Independently on the 
sign of the parameter fi, the leading slope of the 
pulse is steepening due to nonlinearity. In the 
case of the active medium, fi < 0, the ampli- 
tude of the pulse is amplified. In the dissipa- 
tive medium, \i > 0, the amplitude decays. The 
sketch shown on Fig. 4 is valid qualitatively for 
both slow and fast waves. The quantitative dif- 
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ferences arise from different angular dependences 
of the speed of propagation and nonlinear and 
nonadiabatic coefficients (see Fig. 1, 2 and 3). 

6. Astrophysical Applications 

Now, we shall apply the theory developed 
above to two special astrophysical plasma sys- 
tems, a high temperature atomic (T > 10 4 K) 
plasma and the cold (T < 10 3 K) molecular ISM 
gas. 

6.1. High Temperature Gas 

An optically thin plasma in the range of tem- 
perature 10 4 < T < 10 7 K radiates according to 
the relation 

A = p 2 A(^)", (erg cm- 3 s" 1 ) , (60) 

-t i 

where the coefficients A, and the exponent r\ are 
well known in the above range of temperature 
(Vesecky et. al. 1979 ). Therefore, if one as- 
sumes that the plasma heats as ~ p a T b (erg 
cm- 3 s' 1 ) (Dahlburg & Mariska 1988), then, 
it becomes isentropically unstable (p < 0) in 
the ranges of temperature 1.56 x 10 4 < T < 
3.16 x 10 4 K, 2.51 x 10 5 < T < 6.31 x 10 5 
K and 2.0 x 10 6 < T < 3.16 x 10 7 K for 
2.51 xlO 5 <T< 6.31 x 10 5 K for heating by coro- 
nal current dissipation (a = b = 1) and heating 
by Alfven mode/mode conversion (a = b = 7/6), 
and in the range 2.51 x 10 5 < T < 6.31 x 10 5 K 
for a constant heating per unit of mass (a = 1 
and 6 = 0). Thus, in the Sun, the enhance- 
ment of shock waves formation due to thermal 
effects will be relevant in the chromosphere, tran- 
sition region and corona where the corresponding 
e- folding time t e = 2Cg\A\~ 1 , for thermal insta- 
bility is of the order of 5, 2 x 10 2 and 6 x 10 5 s, re- 
spectively. For disturbances with L/Vq « 
the breaking times are of the order of 0.626 1 | 1 
for a magnetic field of about 1 G . This value de- 
creases when the intensity of the magnetic field 
increases. 



On the other hand, Interstellar Medium re- 
gions with temperature in the range 2.51 x 10 5 < 
T < 6.31 x 10 5 K (coronal gas) will be isentrop- 
ically unstable with e-folding times of the order 
of 3 x 10 15 s and a breaking time ti ~ 0.609|/i| _1 
for disturbances with L/Vq ~ VIH> ^ one as " 
sumes a constant heating per unit of mass (a = 1 
and 6 = 0) and a magnetic field Bq of about 
2 x 10~ 6 G. 

6.2. Cold Molecular Gas 

For diffuse and cold (T < 10 3 K) ISM molec- 
ular gas (Van Dishoeck & Black, 1986; Viala, 
1986) one may assume, as a first approximation, 
a heating by cosmic rays and grain photoelec- 
trons (Ibanez and Parravano 1994 and references 
therein), and a cooling (Hollenbach & Mc Kee, 
1979; Hollenbach, 1988) functions to be of the 
form 

T = Hop , (erg cm~ 3 s _1 ) , (61) 

A = 8.5 x W- 5 n 2 L(CO) + n 2 H , 2 L(H 2 ) + 

n H n H2 L(H) , (erg cm' 3 s" 1 ) , (62) 

where n#, n# 2 and n are the number density of 
the atomic, molecular hydrogen, and the total 
number density of particles, respectively. L(CO), 
L(H2) and L(H) are the cooling efficiencies of 
the CO and H 2 molecules and atomic hydrogen, 
respectively. In denser molecular regions other 
heating processes becomes dominant with more 
complex dependence of F on p and T, detailed 
study of which will be reported elsewhere. 

The range of temperature (and particle den- 
sity) at which the gas is isentropically unsta- 
ble (A < 0), is strongly dependent on the ratio 
nn/n and on the energy input parameter Hq. 
In particular, for nn/n = 0.9 these ranges are 
H5 £ T £ 1905 K (136 £ n) and 90 £ T £ 1865 
K (3 £ n) for H = 10~ 2 and H = 10~ 4 , respec- 
tively. However, it is found that the above ther- 
mal instability occurs for high values of nn/n, i.e. 
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for low concentrations of H2 and therefore, when 
the transition from atomic to molecular gas takes 
place. The range of temperature for which the 
gas is thermally unstable decreases when nH 2 / n 
increases and is quenched for high enough values 
of n H Jn (£0.13). 

The ratios between the breaking and the e- 
folding time, %jt e for a disturbance with Gaus- 
sian profile and L/V$ = 10 16 s, have been plotted 
in Figures 5a and 5b, in absence of a magnetic 
field (dash-dot line) as well as for the Galactic 
field value Bq = 2 x 10~ 6 G, for three differ- 
ent values of the angle a, for the fast (a) and the 
slow (b) mode . The breaking times show max- 
ima at a temperature of about 125 K, but the 
breaking times for the magnetoacoustic modes 
are shorter than that corresponding to the hy- 
drodynamic mode. It means that the effect of 
the magnetic field is to decrease the time when 
the breaking of the wave occurs. For the fast 
mode, the effect of increasing the angle between 
the wave number and B is to increase the cor- 
responding breaking time . This dependence re- 
verses for the slow mode (see Figs. 5a and 5b). 

Figure 5c shows the dependence of tb/t e on Hq 
for the fast (dash-line ), slow (continuous line) 
and hydrodynamic (dot-line) mode. The gross 
effect of increasing Hq is to increase the corre- 
sponding breaking times. 

The distance between the point where the dis- 
turbance is originated to the point where the 
shock is formed, is of the order of C s tb. There- 
fore, from above results one may conclude that 
strong inhomogeneities are expected to be formed 
at time scales and scale-lengths shorter than 
those for which the molecular gas becomes ther- 
mally unstable. 

7. Discussion 

We considered dynamics of magnetohydrody- 
namic waves in a nonadiabatic plasma taking 
into account effects of finite wave amplitude. The 
nonadiabaticity was modelled by a general cool- 



ing/heating function of thermodynamic parame- 
ters of the plasma. Depending on a general struc- 
ture of the nonadiabatic function, there are two 
different regimes of the wave evolution in the sys- 
tem, either amplification or decay. The nonlin- 
earity, through generation of highest harmonics, 
leads to the wave steepening and, therefore, to 
formation of shock waves. Our analysis was re- 
stricted by this shock formation time. In the lin- 
early unstable case (usually called the thermal in- 
stability), the formation of the shock takes place 
quicker than in adiabatic case. Consequently, 
nonlinearity depresses the thermal instability of 
magnetoacoustic waves, whose growth is accom- 
panied by the wave steepening and lasts until the 
shock has been formed. 

In the single wave approximation, we have 
derived an evolutionary equation describing de- 
velopment of weakly nonlinear magnetoacoustic 
waves in a weakly nonadiabatic medium. This 
equation is valid for one magnetoacoustic wave 
(either slow or fast) providing that another wave 
is absent from the system. Of course, in real- 
istic situation, both slow and fast waves are ex- 
cited simultaneously. Nevertheless, the difference 
in their propagation speeds gives us a possibility 
to restrict our consideration by the single wave 
approximation, because the waves occupy differ- 
ent locations in time after the excitation (they 
are spatially de-tuned from each other). Since 
the magnetic field brings spatial anisotropy in 
the system considered, parameters of the magne- 
toacoustic wave propagation (the speed, the rate 
of growth or decay, the breaking time) depend 
upon the direction of the wave vector. The third 
magnetoacoustic mode, the Alfven wave, is not 
subject to the influence of the nonadiabaticity 
(as far as the quadratic nonlinear terms approx- 
imation holds). This does not exclude a possi- 
bility of the indirect effect of the nonadiabaticity 
of the plasma on the Alfven wave through non- 
linear coupling with the magnetoacoustic wave 
amplified by thermal instability. The nonlinear 
excitation of the Alfven waves is described by 
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Eq. (|32|). According to above discussion, the ac- 
tive non-adiabaticity of the plasma leads to exci- 
tation of magnetoacoustic waves perturbing the 
variables B x , V x , V z , p and p. Expressing all vari- 
ables through V z by (H)-®), 

we can re-write 

the nonlinear Alfven wave equation ( |32| ) for two 
variables, V y and V z : 



C 2 



v„ 



dV 1 dV l _d_ f dV y 
dz dt dt \ z dz 



2 a 2 (, rn -,av y 



dz 



(63) 

where D t 1 is an inverse differential operator. In 
the case of a harmonical magnetoacoustic wave, 
Eq. d63| ) describes generation of the Alfven wave 
by a parametric instability. But, in the case con- 
sidered, the function V z is prescribed by expres- 
sion (|54|). The spectrum of the exciting magne- 
toacoustic wave is much reacher than harmon- 
ical and evolves in time and space. However, 
qualitatively, we can expect that Alfven waves 
are nonlinearly generated by the thermal insta- 
bility through the discussed mechanism. This 
phenomenon has been omitted from this consid- 
eration. 

The theory developed here may be applied, in 
particular, to an interpretation of origin of den- 
sity, temperature and magnetic field variations 
in astrophysical plasmas, which can be associated 
with the magnetoacoustic weak shock waves (see, 
Krasnobaev, 1975). With respect to the previ- 
ously analysed case of the unmagnetised plasma, 
the presence of the magnetic field gives several 
advantages for the observational registration of 
these phenomena. Firstly, in the vicinity of the 
shock front, the sharp inhomogeneity of the abso- 
lute value of the magnetic field takes place, which 
can be directly registered by analysis of varia- 
tion of parameters of gyro-emission and Zeeman 
splitting. Even if the wave amplitude is suffi- 
ciently small, the sharp gradients in the abso- 
lute value of the magnetic field in the vicinity 
of the shock are accompanied by high densities 



of the electric currents (the transversal compo- 
nent of the current density is proportional to the 
spatial derivative of the magnetic field perturba- 
tion, j y = — (c/^ir)dB x /dz). These current con- 
centrations (or current sheets, in the ideal limit) 
can be regions where particle acceleration takes 
place. Secondly, slow and fast magnetoacoustic 
perturbations, initially excited somehow, propa- 
gate with different speeds and growth rates and, 
therefore, form different spatial structures at dif- 
ferent distance from their source. According to 



formula (48), perturbations of the magnetic field 
are either positively or negatively correlated with 
perturbations of the plasma density, depending 
on the type of the mode (slow or fast), which 
gives us a possibility to distinguish fast and slow 
modes in analysed data. 

From the point of view of the nonlinear wave 
dynamics, propagation of the MHD waves in the 
optically thin plasma is a very attractive prob- 
lem, because the plasma is acting in this case as 
an active medium, supplying energy to the waves. 
The simplest model considered here can be fur- 
ther developed. For example, high frequency dis- 
sipation (modelled by introduction of a Burgers 
term in equation (|52|)), will compete with the 
amplification and depress nonlinear generation 
of highest harmonics. When these phenomena 
are in balance, there exist dissipative stationary 
waves with "saw-teeth" shape. Also, the gen- 
eration of the second harmonics and so the en- 
ergy transfer to small scales can be depressed 
by nonlinear dispersion (Whitham, 1974). This 
process is based upon either cubic nonlinear- 
ity, or quadratic nonlinearity through cascade 
self-interaction. The last mechanism has two 
stages: nonlinear generation of second harmon- 
ics and interaction of the second harmonics with 
a backward wave. In the case of magnetoacous- 
tic waves, both mechanisms are acting together 
(Nakariakov et al., 1997b) leading to the appear- 
ance of a cubically nonlinear term in the evolu- 
tionary equation. Obviously, the cubically non- 
linear effects should be taken into account only 
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when quadratically nonlinear processes are sup- 
pressed due to, e.g., linear high frequency disper- 
sion. Effects of dispersion (connected possibly 
with electron inertia or, which is more impor- 
tant for astrophysical applications, with inhomo- 
geneity of the plasma) can be modelled by intro- 
duction of a Korteweg - de Vries (or Benjamin 
- Ono or Leibovich - Roberts, depending on the 
specific situation) term in evolutionary equation 



(52). Particularly, the dispersion of the medium 
can give a possibility for existence of autosolitons 
and stochastic dynamical regimes (see, e.g. En- 
gelbrecht, 1989, Nakariakov and Roberts 1999). 
These phenomena can be registered in astrophys- 
ical plasmas and be applied to interpretation of 
the observations. 

Finally, we must remark that in atomic as well 
as in molecular plasma occurring in astrophysics, 
many other physical effects likely go into play, 
such as inhomogeneity in the background due 
to self-gravitation and/or to radiation transport, 
and dissipative effects ( viscosity, thermal con- 
duction, ionization-recombinations and chemical 
reactions). These phenomena have to be taken 
into account in more developed models. How- 
ever, the above results indicate that inside ther- 
mally unstable plasma regions, substructures can 
be originated by nonlinear effects in times shorter 
than an e-folding time for evolving the isentropic 
thermal instability. 



8. Acknowledgements 

VMN is grateful to Prof Bernard Roberts and 
Dr Mikhail Ruderman for a number of stimu- 
lating discussions. CAMB and MHIS thank to 
the CDCHT of Universidad de los Andes for its 
support through projects C-931-99-05-B and C- 
912-98-B-05, respectively. 



REFERENCES 

Balbus, SA. and Soker, N., 1989, ApJ, 341, 611. 

Bodo, G., Ferrari, A., Massaglia, S. and 
Tsinganos, K., 1985, A&A, 149, 246. 

Dahlburg, R.B., De Vore, C.R., Picone, J.M., 
Mariska J.T. and Karpen, J.T., 1987, ApJ, 
315, 385. 

Dahlburg, R.B. and Mariska, J.T., 1988, Solar 
Phys., 117, 51. 

Engelbrecht, J. (Ed.), 1989, Nonlinear Waves in 
Active Media, Springer- Verlag, London. 

Field, G.B., 1965, ApJ, 142,531. 

Flannery, B.P. and Press, W.H., 1979, ApJ, 231, 



Heyvaerts, J., 1974, A&A, 37, 65. 

Hollenbach, D. & Mc Kee, C. 1979, ApJS, 41, 
555. 

Hollenbach, D. 1988, Astrophys. Lett. Commun., 
26, 191. 

Ibahez S., M.H., 1985, ApJ, 290, 33. 

Ibahez S., M.H. and Escalona T., O.B., 1993, 
ApJ, 415, 335. 

Ibahez S., M.H. and Parravano, A., 1983, ApJ, 
275, 181. 

Ibahez S., M.H. and Parravano, A., 1994, ApJ, 
424, 763. 

Kalkofen, W., 1997, ApJ. Lett, 486, L145. 

Karpen, J.T., Antiochos, S.K., Picone, J.M. and 
Dahlburg, R.B., 1989, ApJ, 338, 493. 

Krasnobaev , K.V.1975, Izv. Akad. Nauk SSSR, 
Mech. Zhidk. Gaza 2, 126. 

Krasnobaev, K.V. and Tarev, V.Y., 1987, Sov. 
Astron. 31, 635, (KT). 



12 



Kruskal, M. and Schwarzschild, M.,1966, in Mag- 
netohydrodynamics Stability and Thermonu- 
clear Containment, Eds. A. Jeffrey and T. 
Taniuti, p. 63 (Academic Press, New York) 

Landau, L.D. and Lifshitz, E.M.,1987, Fluid Me- 
chanics (Pergamon, London). 

Nakariakov, V.M. and Oraevsky, V.N., 1995, So- 
lar Phys., 160, 289. 

Nakariakov, V.M. and Roberts, B., 1999, Physics 
Letters A 254, 314. 

Nakariakov, V.M., Roberts, B. and Murawski, 
K., 1997a, Solar Phys. 175, 93. 

Nakariakov, V.M., Roberts, B. and Petrukhin, 
N.S., 1997b, J. Plasma Phys. 58, 315. 

Nakariakov, V.M., Roberts, B. and K. Murawski, 
K., 1998, A&A, 332, 975. 

Oppenheimer, M., 1977, ApJ, 211, 400. 

Parker, E.N., 1953, ApJ, 117, 431. 

Priest, E.R., 1982, Solar Magneto Hydrodynam- 
ics (Reidel, Dordrecht), pp. 235, 278, 382. 

Reale F., Serio, S. and Peres, G., 1994, ApJ, 433, 
811. 

Rosner, R., Tucker, W.H. and Vaiana, G.S.,1978, 
ApJ, 220, 643. 

Tarev, V.Y., 1993, Astron. Lett. 19, 441. 

Van Dishoeck, E. F. and Black, J. M., 1986, 
ApJSS, 62, 109. 

Viala, Y. P. 1986, A&ASS, 64, 391. 

Vesecky, J.F., Antiochos, S.K. and Underwood, 
J.H., 1979, ApJ,233, 987. 

Weymann, R., 1960, ApJ, 132,452. 

Whitham, G.B., 1974, Linear and Nonlinear 
Waves (Wiley- Interscience Pub., New York) 

This 2-column preprint was prepared with the AAS IATeX macros v4.0. 



13 



Fig. 1. — The linear speeds C/C s as a func- 
tion of a for different values of (3. Continuous 
and dashed lines correspond to the fast and slow 
modes, respectively. 

Fig. 2.— (a) The non-adiabaticity (^C 2 S \A\- V ) 
and (b) the nonlinearity e parameters as func- 
tions of a for different values of (3 corresponding 
to the fast wave. 

Fig. 3. — As Figure 2 for the slow wave. 

Fig. 4. — Development of the shape of the initial 
Gaussian magnetoacoustic pulse with the width 
L situated initially t = in z = 0, and the indi- 
cated time in units of L/Vq. In the active medium 
(a) (/i < 0) the pulse is amplified. In the dissi- 
pative medium (b) fi > the pulse decays. 

Fig. 5. — The breaking times for the fast (a) and 
slow mode (b) as a function of temperature T, 
for H = 1(T 3 , B = (dash-dot line), B = B G , 
ntj/n = 0.9 and three values of the angle a = 
it/36 (dash line ), 7r/4 (solid line) and 177r/36 ( 
dot line). Fig.5c is as Fig. 5a for two different 
values of H , 1(T 3 and 1(T 4 . 
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